Stability of symmetric tops via one variable calculus 
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^ Abstract 

00 ! 

^SJ , We study the stability of symmetric trajectories of a particle on the Lie group S0{3) whose 

motion is governed by an 50(3) x 50(2) invariant metric and an 50(2) x 50(2) invariant 
potential. Our method is to reduce the number of degrees of freedom at singular values of 
the 50(2) X 50(2) momentum map and study the stability of the equilibria of the reduced 
systems as a function of spin. The result is an elementary analysis of the fast/slow transition in 
' the Lagrange and Kirchhoff tops. More generally, since an 50(2) x 50(2) invariant potential 

00 . on 50(3) can be thought of as Z2 invariant function on a circle, we get a condition on the 

' second and fourth derivatives of the potential at the symmetric points that guarantees that the 

, corresponding system gains stability as the spin increases. 

1 Introduction 

■ The motivation for writing this paper is two- fold. The first one is to see the stable/unstable 

transition in the Lagrange top without resorting to the machinery of equivariant singularity theory 
as m 0, jCvdMl let alone the more sophisticated techniques of fLRSMH . The idea is to exploit the 
\^ • SO{2) X 5*0(2) symmetry and use reduction instead. The resulting system can then be analyzed 



by looking at the signs of two derivatives of a function of one variable at one particular point. The 
second motivation is to explain why it is that both the Lagrange top and the Kirchhoff top [BZ] 
undergo the Hamiltonian Hopf bifurcation as their spin varies. 

I will show that by carrying out reduction at certain singular values of the S0{2) x S0{2) 
momentum map one can understand the bifurcation in both systems (in fact in an infinite family 
of systems) as a phenomenon of the following sort. Consider a particle on a line whose behavior is 
governed by a 1-parameter family of potentials 

/a(x) = ^x^-x2 + aV, 
o 

where A is the parameter. Then for |A| < 1 the origin is a local maximum of the potential, hence 
an unstable equilibrium of the system. Global minima occur at x = ±(1 — A^)^/^. As |A| increases, 
these stable points move closer to the origin until they collide at |A| = 1; the origin becomes a 
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a global minimum (hence a stable equilibrium). In other words, the reduced system undergoes a 
generic Hamiltonian Z2 invariant bifurcation — the figure eight bifurcation (cf |GoSt] ] ) . 

There is another interpretation of the results of this paper. Namely, the systems that are 
being studied are three degrees of freedom completely integrable systems. Thus the nondegenerate 
singularities of the corresponding Lagrangian foliation are, in Eliasson's classification of elliptic 
or of hyperbolic type. It turns out that both types of singularities are present in our systems 
and that the closures of the sets of elliptic and hyperbolic points intersect at some degenerate 
singular point. This fact was already observed by Eliasson in the case of the Lagrange top (op. 
cit.). The point I would like to make is that the Hamiltonian Hopf bifurcation in these systems 
can be described as a transition from an elliptic to a hyperbolic singularity through a degenerate 
singularity. It would be interesting to see if such transitions are in some sense generic for the 
three degrees of freedom completely integrable systems. In higher dimensions many more types of 
nondegenerate singularities of the Lagrangian foliation can occur. This suggests that transitions 
between these singularities should give rise to interesting bifurcations. 

2 The set-up and the first 50(2) reduction 

The configuration space of the systems we will study is the Lie group S0{3) and the Hamiltonians 
H are of the form Kinetic Energy + Potential Energy, where the kinetic term comes from a left 
SO (3) and right S0{2) invariant metric, i.e. the moments of inertia satisfy Ii = I2 ^ I3, and the 
potential is left-right 50(2) invariant. Here 50(2) is identified with the isotropy group of 63, the 
third vector in the standard basis {ei, 62, 63} of R'^. For example l^(^) = (^63, 63) is the potential 
of the Lagrange top and ^^(A) = {Ae^, ei)^ + {Ae^, 62)^ + c(Ae3, 63)^, c a positive constant, is the 
potential of the Kirchhoff top. 

Since the left action of the group 50(2) on the cotangent bundle of 50(3) is free, we may carry 
out the Marsden-Weinstein-Meyer reduction. It amounts to fixing the rate of rotation of the top 
about its axis of symmetry. A routine calculation shows the reduced spaces, as manifolds, are the 
cotangent bundle of the two sphere 5^ with the reduced symplectic structures of the form 

ojx = uj'pt 52 + Xu, 

where A is the value of the 50(2) momentum map, ujrp*s^ is the standard symplectic form on the 
cotangent bundle of the sphere and v is the standard area form on the sphere (pulled up to the 
cotangent bundle). For example, to compute these reduced spaces one could use Kummer's results 
[K] or the fact that 50(3) and the Euclidean group £'(3) form a dual pair. The form ly is often 
referred to as the magnetic form. This calculation will not be carried out in the paper. The reader 
may simply assume that we are studying a family of Hamiltonian systems on (T*S'^,uj\), A G R. 

There is one action of 50(2) left on the reduced space. If we identify the sphere with the 
standard round sphere in R'^ then the action under the identifications made above is the lift of 
the rotation about the 63 axis. The reduced Hamiltonians are again of the form Kinetic Energy + 
Potential Energy, where the kinetic term comes from a round metric on the sphere and the potential 
V is an 50(2) invariant function. For example the Lagrange top potential is V{x) = X3\g2 and the 
Kirchhoff top potential is V{x) = (x^ + + C3;3)|5'2 = (1 + (c — 1)3:3)152. Here (xi,X2,X3) are the 
standard coordinates on R'^. 

The momentum map ^\ : {T* S'^ ,uJx) — > R is the sum of two momentum maps: one coming 
from the lifted action of 50(2) on {T* S"^ , uj^* s^) ^'^'^ the other coming from the action of 50(2) 
on (5^, Az^). We fix the constant in the definition of the momentum maps by requiring the map to 
be at the North pole. 
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We think of the reduced systems {T*S'^,u}\, h) as two degrees of freedom systems depending on 
one real parameter A, the spin. The action of 50(2) fixes the North and South poles. Consequently 
these points are critical for our Hamiltonians h. We would like to understand how the behavior of 
the systems near the symmetric points changes as we vary the parameter A. We will consider the 
North pole. The arguments regarding the South pole are exactly the same. In fact the Hamiltonian 
for the Kirchhoff top is invariant under the Z2 symmetry interchanging the poles. 

To understand the behavior of the system near the North pole we reduce with respect to the 
remaining symmetry at ^\ = 0, the value of the normalized momentum map at the North pole. 
Since the North pole is a fixed point, the standard Marsden-Weinstein-Meyer procedure does not 
apply. We carry out the reduction in the next section. 

3 Reduction at a critical value of the momentum map 

Let us briefly summarize the procedure that we are going to follow. A full account is given in 
HSL| . Our point of view is that given a Hamiltonian action of a compact Lie group K on a. 
symplectic manifold (M, w) with momentum map F : M — > k*, the reduced space at a coadjoint 
orbit Oa = -fC • a C k* is the Hausdorff topological space 

■.= F-\0^)/K. 

Smooth K invariant functions on M descend to continuous functions on M^; we denote these 
functions by C°°(Mq). Algebraically 

C°°(M«) = C°^(M)^/I„, 

where C°°{M)^ denotes the algebra of smooth K invariant functions and /q is the ideal in C°°(M)^ 
consisting of functions that vanish on F~^{Oa)- It is not hard to see that la is a Poisson ideal 
[ ACG[| . Indeed, the Hamiltonian flows of K invariant functions preserve the set F-i(0 a), SO for 
any / G C°°{M)^ and any h ^ la we have {f,h} = 0. Therefore C°°(Mq,) is a Poisson algebra. 
The Poisson bracket allows us to associate flows to functions. It also makes it possible for us to say 
when two reduced spaces are isomorphic. 

Definition 3.1 Two reduced spaces X and Y are isomorphic if there exists a homeomorphism 
(J) : X ^ Y such that the induced map (j)* : C°°{Y) — > C°°(X) is a Poisson isomorphism. 

The following example will be very important for us. In fact, thanks to the equivariant Darboux 
theorem, this is precisely the situation that we are interested in. Consider the standard lifted action 
of S0{2) on the cotangent bundle r*R^. Let {x,y) be the coordinates on and {x,y,px,Py) the 
canonical coordinates on r*R^. The standard symplectic form uj is given in these coordinates by 
CO = dx A dpx + dy A dpy and the momentum map F by 

F{x,y,px,Py) = xpy - ypx- 

The map tf^ : T*R T*'R?, {u,pu) ^ (n, 0,pu,0) is a symplectic embedding whose image lies 
entirely in the zero level set F-i(O) of F. Moreover, all the orbits of S0{2) in F-i(O) (except {0}) 
intersect the image in exactly two points: (u, 0,p„, 0) and {—u, 0, —pu, 0). Consequently tp descends 
to a homeomorphism (j) : T*R/Z2 — > F^^{0)/ S0{2) = (T*R^)o where Z2 acts on r*R by the map 

{u,Pu) ^ (-M, -Pu)- 

Note that since any 5*0(2) invariant function on T*R^ pulls back via V to a Z2 invariant 
function on T*R, (p pulls back smooth functions on the reduced space {T*'R?)q to smooth functions 
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on the orbifold T*R/Z2. Since ip is a symplectic embedding the puh-back (p* is Poisson. It is easy 
to see that (j)* is injective: i;^ is a homeomorphism. To prove surjectivity of the puh-back we use a 
theorem of G. Schwarz on invariant smooth functions [Sch]. The Z2 invariant polynomials on T*R 
are generated by three polynomials: u'^, and upu- Schwarz's theorem in this case says that any 
smooth Z12 invariant function is a composition of a smooth function on R'^ with the map T*R R^, 
{u,Pu) ^ iu^,pl,upu), i.e., if / G C°^{T*'R)^^ = C°^(T*R/Z2) then there is / G C°^(R3) with 

f{u,Pu) = f{u^,pl,upu). 

Now = il)*{x^ + y^), p^ = ip*{p'^ ~^ Py) ^"^^ "^P^ ~ '^*i^Px + UPy)- Therefore any function 
/ G C°°(r*R/Z2) is in the image of (/>*. 

Conclusion: the reduced space (r*R^)o is isomorphic to the orbifold T*R/Z2 (the cotangent 
bundle of a line reduced by a finite group Z2). 



3.1 Stability 



Marsden and Weinstein | MW ] proved that stable equilibria of reduced systems correspond to stable 
relative equilibria of the original system. More precisely they considered a Lie group G acting on 
a symplectic manifold (M, w) in a Hamiltonian fashion with corresponding momentum map $ : 
M — > g* and an invariant Hamiltonian h G C°°(M)'^. For a regular value a of the momentum map 
let ha denote the corresponding reduced Hamiltonian on the reduced space Ma ■= ^~^{G • a)/G. 
Marsden and Weinstein showed that if x G Ma is a critical point of ha and the Hessian of ha at 
X is definite then the trajectory of the original Hamiltonian h through a point m G <l>~^(a) that 
corresponds to x G Ma is stable relative to the group action. In fact the result holds under a 
slightly weaker hypothesis, namely that the action of the group G is free at m, i.e., the reduced 
space Ma is smooth near x. Thus a being a regular value is not really needed. 

If a point m is symmetric, i.e., has a nontrivial isotropy, then it is automatically a critical point 
of any invariant Hamiltonian. Since it is also a critical point of the momentum map Marsden and 
Weinstein's result does not apply. However, for the action of SO{2) on the cotangent bundle of R^ 
that we are considering, the correspondence still holds. 

Lemma 3.2 Suppose a function h G C°°(T*R^) is S0{2) invariant. Then if the Hessian of 
ho := ^|{(u,o,pu,o)gt*r2} is definite at the origin, the origin = (0,0,0,0) G T*'R? is S0{2) stable, 
i.e., given any S0{2) invariant neighborhood U of in the cotangent bundle T*R^ there is an 
SO{2) invariant neighborhood VofO such that for all x & V the trajectory of the Hamiltonian flow 
of h through x stays in U. 

Proof. It is no loss of generality to assume that the Hessian 'H(/io) of /iq at is of the form 

H{ho) = au^ + bp'^; a,b > 0. 
Indeed, any quadratic form on R^ can be diagonalized by an orthogonal change of coordinates, 
1^1 c s \ I u \ where + = 1. Moreover, this change of coordinates is induced 

Pu J \ -s c J y pu J 

by an S0{2) equivariant change of coordinates on T*R^; it is given by the matrix 



c 

-s 




s \ 

c s 

c 

-s c J 
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Since the Hessian of ^ at is SO{2) invariant and since o,p„,o)Gr*R2} = 'H{ho), the Hessian 

of h has to be of the form 

n{h) = a{x^ + y2) + b{pl +pI) + cF, 

where c is a constant and F = xpy — ypx is the momentum map. It cannot have any more terms 
for then the Hessian of the reduced Hamiltonian H{ho) would have more terms. 

Let p be a a compactly supported smooth function on the cotangent bundle of which is 
identically 1 near 0. Then near the origin the Poisson bracket {h, pcF} is zero. The Hessian of 
h := h — pcF is positive definite at 0, hence the origin is a stable equilibrium of h. On the other 
hand, since the flows of h and of pcF commute near the origin, the flow of the original Hamiltonian 
is the composition of the flows of h and of cF. Therefore G r*R^ is SO{2) stable for the flow of 
h. □ 



3.2 Reduction of the cotangent bundle of the 2-sphere 

Consider coordinates on the upper hemisphere of 5^ induced by the projection {x,y,z) ^ {x,y). 
The inverse map (p is given by 4>{x,y) =i— > (x, y,\/l — x"^ — y^). These coordinates allow us to 
compute the most interesting part of the reduced system. We don't need the coordinates to com- 
pute the reduced space per se, rather the coordinates are convenient for computing the reduced 
Hamiltonians. Since ^ 

(p* {xdy hdz + ydz Adx + zdx A dy) = -dx A dy, 

where z = \fY^^x^^'^^y^ , the symplectic form bj\ in the canonical coordinates corresponding to 
(x, y) is 

u\ = dx A dpx + dy A dpy H — dx A dy. 
The momentum map $a in these coordinates is then 

*A(a:,y,Pz,]3j/) = yPx - xpy - Xz + X. 

Note that the North pole in these coordinates is the origin and that 'I>(0) = 0. 

The metric g by assumption is the one that comes from the embedding of S'^ as a round sphere 
in R^. We now normalize the metric by setting the radius of the sphere to 1. Then the sphere is 
cut out by the equation 

x'^ + y'^ + z'^ = 1, (1) 

In coordinates the metric g is given by 





Consequently 

We are now in position to compute the reduced Hamiltonian. From the discussion in the 
beginning of the section we know that locally the reduced space is the cotangent bundle of the 
interval (—1,1) divided by the involution {u,Pu) ^ {—u,—pu)- Since our coordinates are not 
Darboux we proceed as follows. Set y = 0. Then ^\{x, 0,Px,Py) = if and only if 

Py = A(l - z)/x. 
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Note that for y = 0, 

l-z = l- =l-{l-l-x^- + 0{x^)) = Ix^ + Ix^ + 0{x^), 

2 8 2 8 

so 

Py{x) = -{x + -x^ + 0{x^)) 

is smooth at x = 0. Consider now a map ip : T*(— 1,1) —>■ $^''^(0), {u,Pu) >—>■ (u, 0,Pu,A(l — z)/u) 
where now z = \/\ — v?. As in the example considered above, the composition 

r*(-i,i)^$^i(o)^$^i(o)/so(2) 

is a branched double cover. Therefore, we will simply work on r*(— 1, 1) keeping the involution in 
mind. Now the pull-back by '0 of the kinetic energy part of the Hamiltonian is 

where G{u) = (1 — u^)- The potential energy part pulls back to an even function of «, hence to a 
function of u^. For example, for the Lagrange top the corresponding function is 

V{u^) = ^/l-v? 

and for the Kirchhoff top 

V{u^) = 1 + (c- 1)(1 - n^). 
The critical points of the reduced Hamiltonian 

are in one-to-one correspondence with the critical points of the effective potential 

Moreover, relative maxima of Ux correspond to unstable equilibria of Hx and relative minima to 
stable equilibria. 

Remark Since we are working on the branched double cover of the reduced system, any pair of 
critical points of Ux of the form ±u, u ^ correspond to the same critical point of the reduced 
system. Therefore only the non-negative critical points of the potential need be considered. 

3.3 Critical values of the effective potential Ux{u), u>Q 
Since 

there exists a change of variables u = t{v) so that 

t{v) ) ■ 
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Note that t{v) is an odd function of the form 



t{v) = 2v- 2v^ + 0{v^). 

Then the transformed reduced potential V{t{v)'^) is even in v and so has to be of the form 
V{t{v)'^) = /(f ^) for some smooth function /. Therefore the study of the bifurcation of the critical 
points of the potential U\{u) for u small and non-negative is reduced to studying the bifurcation 
of the critical points of Fx{v) = ^v'^ + /(v^) for v small and nonnegative. 

It will be useful to express the derivatives /'(O) and /"(O) in terms of the derivatives V'(0) and 
V"{0). Since t{v)^ = Av^ - Sv"^ + 0{v^) and T(t;)^ = IGv"^ + 0{v'^), we have 

V{T{vf) = V{0) + V'{0)T{vf + V"{0)r{v)'^ + 0{v^) 

= V{0) + 4V'{0y + 8(F"(0) - V'{0)y + 0{v^). 



Therefore 



/'(O) = AV'iO) and /"(O) = 8(F"(0) - y'(O)). (2) 



3.4 Bifurcation of critical points of Fx{v) — + /(w^), v>^ 

It is no loss of generality to assume that A > 0. We also make the genricity assumptions: /'(O) ^ 
and /"(O) / 0. 

Thus if /'(O) > then w = is the only critical point for all values of A, i.e., no bifurcations occur. 
Physically this mean that if a top in a straight up position with zero spin is stable, then it is stable 
in this position for all values of spin. 

Assume now /'(O) < 0. Then /'(O) = -c? jl for some a > 0. Consider A) = A^ + 2/'(t). 
Since {f'y{0) 7^ 0, by the implicit function theorem there exists a smooth function g{\) defined 
near A = a such that g{a) = and 4>{g{X),X) = 0, i.e., 

^A^ + /'(5(A)) = 0. (3) 



Note that since 



we have 



d 



X=a 



</.(5(A),A)=[A + r(5(A))5'(A)]|,^„ 



f'Wia) = -a < 0. 

Thus if /"(O) > then g{X) is strictly decreasing near A = a. Consequently there are e, (5 > so 
that for < V < 5 the equation A^ + /'(f ^) = has no solution whenever a < A < a + e and for 
a — e < X < a the equation has exactly one solution, namely {v, A) = (51(A) A); see figure 1. 

Similarly if /"(O) < then g{X) is strictly increasing near X = a. Consequently there are e,S > 
so that ioT < V < S the equation A^ + f'{v^) = has no solution whenever a — e < X < a and for 
o < A < a + e the equation has exactly one solution, namely (f , A) = (5(A) A). 

Finally we determine the stability of these critical points. Since 

l^Fx{v)=X^ + 2r{v^) + Av^r{v^) 
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/"(O) > /"(O) < 
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Figure 1: bifurcation of critical points of Fx{v) 



we have, using (3), 



and 



Fx{v)=4\g{\)\f"{g{\)) which is 



^=5(A)i/2 



> if /"(O) > 
< if /"(O) < 



Fx{v) = 2(A2 + /'(O)) which is 



v=0 



> if A > v /-2f(0) 
< if A < V-2/'(0) 



Conclusions Assume that /'(O) < 0. 

If /"(O) > then as A decreases below Aq = v^— 2/'(0), f = changes from a local minimum to a 
local maximum and a one parameter family of local minima v = ^(A)^/^ bifurcate from v = 0. 

If /"(O) < then the situation is reversed. For A < Aq and v small, the origin is a unique 
critical point, and it is a local maximum. As A increases past the critical value, v = becomes a 
local minimum and a one parameter family v = g{\)^^^ of local maxima bifurcate from it. 



8 



These results easily translate into a description of the bifurcation of nonnegative critical points 
of the potential 

Ux{u) = -[^ n ) +^(-)' 

where without loss of generality we shall assume that A > 0. Recall that 

/'(O) = 4V'{0) and /"(O) = 8{V"{0) - V'{0)). 

It follows that if V^'(O) > then u = is a local minimum for all values of A and no bifurcations 
occur. We now assume that V'{0) < 0. Then there are two alternatives. 



1. If V"{0) > V'{0) then as A decreases below Aq = \/—V'{0), u = changes from a local minimum 
to a local maximum of the potential Ux and a 1-parameter family u = p{X) of local minima 
bifurcate off from u = 0. Here p{X) is smooth for A < Aq and Holder 1/2 at Aq- Consequently the 
corresponding system on the cotangent bundle of the two-sphere undergoes a Hamiltonian Hopf 
bifurcation. 

2. If V"{0) < V'{0) then for < A < Aq and u small, the point n = is a unique local maximum. 
As A increases past Xq, u = becomes a local minimum and a 1-parameter family u = p{X) of local 
maxima bifurcate from it. 

For the Lagrange top y(n2) = ^/l — V?, so the first alternative holds with Aq = 2 (this agrees with 
1^). For the Kirchhoff top V{v?) = 1 + (c — 1)(1 — v?). So if c < 1 there is no bifurcation and if 
c > 1 the first alternative holds. 
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